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ABSTRACT 

The intermediate- to high-mass star-forming region IRAS 20343-1-4129 is an excehent 
laboratory to study the influence of high- and intermediate-mass young stellar ob- 
jects on nearby starless dense cores, and investigate for possible implications in the 
clustered star formation process. We present 3 mm observations of continuum and 
rotational transitions of several molecular species (C2H, C-C3H2, N2H+, NH2D) ob- 
tained with the Combined Array for Research in Millimetre-wave Astronomy, as well 
as 1.3 cm continuum and NH3 observations carried out with the Very Large Array, 
to reveal the properties of the dense gas. We confirm undoubtedly previous claims of 
an expanding cavity created by an ultracompact Hii region associated with a young 
B2 zero-age main sequence (ZAMS) star. The dense gas surrounding the cavity is 
distributed in a filament that seems squeezed in between the cavity and a collimated 
outflow associated with an intermediate-mass protostar. We have identified 5 mil- 
limeter continuum condensations in the filament. All of them show column densities 
consistent with potentially being the birthplace of intermediate- to high-mass objects. 
These cores appear different from those observed in low-mass clustered environments 
in sereval observational aspects (kinematics, temperature, chemical gradients), indi- 
cating a strong influence of the most massive and evolved members of the protocluster. 
We suggest a possible scenario in which the B2 ZAMS star driving the cavity has com- 
pressed the surrounding gas, perturbed its properties and induced the star formation 
in its immediate surroundings. 
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1 INTRODUCTION 

Most of the stars of all masses in the Galaxy form in rich 
clusters. Despite this, the details of the clustered star forma- 
tion process are still poorly understood. Studies of low-mass 
protoclusters have started to unveil similarities and differ- 
ences between isolated a nd clu stered dense coresfe.g. Andre 
et al. l2007l , Foster et al. |2009| . Friesen et al. 120091 '). Globally 
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these studies suggest that cluster environment has a rela- 
tively smaller influence on the properties of the cores (tem- 
perature, mass, velocity dispersion, chemical abundances of 
early phase molecules) than is typically assumed (Foster et 
al. 2009). However, the conclusions described above do not 
include observations of high-mass star forming regions. Be- 
cause the phenomena associated with massive star formation 
have a stronger impact on the environment (massive out- 
flows, UV radiation, expanding Hll regions), it is plausible 
that these energetic phenomena have major effects on the 
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surrounding dense material. The study of such interaction 
is especially important to quantify the effect of protostellar 
feedback on the environment and test recent models of high- 
mass star formation including outflows and radiation from 
the n ewly b orn stars (e.g. Krumholz et al. 2011, Hennebelle 
et al. l201lD . 

Interferometric observations of dense gas and dust trac- 
ers (N2H^, mm continuum, NH3) have revealed the presence 
of pre-stellar core candidates surrounding ultracompact Hll 
regions (UC Hll's) and other massive young stellar objects 
(YSOs) that do show evidences of such an interaction. For 
example, Fontani et al. (2009) found that in the protoclus- 
ter associated with IRAS 05345+3157 the kinematics of two 
pre-stellar core candidates is influenced by the passage of a 
massive outflow. UV radiation and powerful outflows affect 
the chemistry of starless cores in IRAS 20293+3952 (Palau 
et al. 2007a). On the other hand, a crucial chemical process 
in pre-stellar cores, i.e. the deuteration of species like N2H"'" 
and NH3, seems to remain as high as in pre-stellar cores 
isolated and a ssocia ted with low-mass star forming region s 
(Fontani et al. |2008| . Busquet et al. 2010, Pillai et al. bOllT ). 
Therefore, to date it is not clear if and how the presence of 
massive objects affects the properties and evolution of the 
other (pre-)protocluster members. 

The protocluster associated with IRAS 20343+4129 
(hereafter 120343) represents an excellent laboratory to 
study this issue. The IRAS source is located to the north- 
eastern side of the Cygnus OB2 association, at 1.4 kpcofdis- 
tance from the Sun (Sridharan et al. 2002, Rygl et al.|2oi3), 
and two bright nebulous stars, IRS 1 (north) and IRS 3 
(south), are found inside the IRAS error ellips e whe n ob- 
served at high angular resolution (Kumar et al. 120021 ). The 
bright infrared stars are embedded in a cometary-like cloud 
whose head, facing the Cygnus OB2 association, is bright at 
centimeter wavelengths and whose tail, bright in the mid- 
infrared, is extending for about 10' (~ 4 pc) towards the 
north-east (Fig. [1} . This kind of clouds are also known as 
bright rimmed clouds. 

Thanks to interferometric observ ations of ^^CO and 
1.3 mm continuum, Palau et al. (J2007bj) concluded that IRS 
1 is an intermediate-mass Class I YSO driving a molecular 
outflow in the east-west direction, while IRS 3 is likely a 
more evolved intermediate/high-mass star. This is further 
confirmed through mid-infrared photometric an d spe ctro- 
scopic observations, from which Campbell et al. (|200a ) also 
estimated a bolometric luminosity of the order of 1000 Lq 
for both IRS 1 and IRS 3. Furthermore, IRS 3 is at the cen- 
tre of an UC Hll-region detected through VLA centimeter 
continuum emission (Carral et al.[l999j), and of a fan-shaped 
emission in the 2.12 um_rovibrational line of molecular hy- 
drogen (Kumar et al .l2002l) . East and west of this fan-shaped 
feature, Palau et al. l|2007a l detected molecular gas and dust 
resolved int o severa l millimeter continuum compact sources. 
Palau et al. l|2007br ) interpreted these starless condensations 
as being accumulated on the walls of the expanding shock 
front, but could not derive firm conclusions on their origin 
and nature. 

This work aims at better understading the nature of 
the dense cores in 120343, and its relation with the neigh- 
bouring more evolved objects. To achieve the goal we per- 
formed observations of molecular species obtained at high 
angular resolution with the Combined Array for Research 



in Millimeter Astronomy (CARMA) at 3 mm and the Very 
Large Array (VLA) at 1.3 cm. All selected molecular transi- 
tions are commonly used to characterise dense gas: (i) NII3 
and N2H^ are excellent tracers of dense and cold cores be- 
cause either species do not suffer from depletion up to ~ 10^ 
cm"'', and NH3 is extensively used as thermometer in both 
low- and high-mass star forming regions; (ii) NII2D provides 
an estimate of the degree of deuteration (with NH3). This 
combination of diagnostic lines was successfully used by Bus- 
quet et al. (2010) to identifying pre-protostellar cores in the 
proto-cluster associated with IRAS 20293+3952; (iii) C2H 
and C-C3H2 are both high-density PDR tracers useful to 
shed light on the interaction among the cold gas and the 
UV radiation field coming from IRS 1 and IRS 3. C2II is 
also a tracer of cold gas (e.g. Beuther et al. 2008, Padovani 
et al. '2009l). In this paper we concentrate on the gas mor- 
phology, temperature and kinematics of the region adjacent 
to IRS 1 and I RS 3, and confirm the hypothesis proposed 
by Palau et al. (|2007bj) that IRS 3 is opening a cavity in the 
surrounding dense gas and starless material is being accu- 
mulated on the cavity walls. In Sect.[5]we describe the obser- 
vations. The observational results are presented in Sect. |3l 
and discussed in Sect, [l] In Sect. [5] we summarise the main 
findings of the work and give a general conclusion. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 CARMA 

3 mm CARMA observations of 120343 were obtained on 29 
Mar 2010 in C- and 01 May 2010 in D-configuration under 
good weather conditions for observations at 3 mm, charac- 
terized by about 5 mm of precipitable water and atmospheric 
noise rms of about 300 /im as measured on a baseline of 
100 m at the frequency of 2 25 GHz . The phase centre was 
the same as in Palau et al. (|2007bl ), namely: R.A. (J2000) 
= 20''36'"0r:3 and Dec. (J2000) = 41°39'57."2. The local 
standard of rest velocity of the cloud is assumed to be 11.5 
km sT^ , as determined from single-dish ammonia observa- 
tions (Sridharan et al. 2002). The primary beam of the 10 m 
and 6 m dishes at about 85 GHz is ~ 73 " and 121 ", respec- 
tively. The single-side-band system temperature during the 
observations was below 150 K. During C-configuration ob- 
servations, the correlator provided 4 bands which were con- 
figured to observe the continuum, the C2H, the NH2D and 
the N2H"'" lines simultaneously. D-array observations were 
obtained with the new CARMA correlator, which provides 
more bands. Two 500 MHz bands were used to observe the 
continuum and 5 bands set up to observe the C2H, N2H+, 
NH2D, CCS and C-C3H2 line emission. The pass-band was 
calibrated by observing 1733-130; flux calibration was set 
by observing MWC349. The estimated uncertainty of the 
absolute flux calibration is 10% , and it is determined from 
periodic observations of MWC349. Atmospheric and instru- 
mental effects were corrected by observing the nearby quasar 
2007+404 every 15 minutes. 

The tracers observed and the main observational pa- 
rameters (frequency, synthesised beam, linear resolution, 
spectral resolution, la rms channel noise, largest detectable 
angular scale) are reported in Table [1] The CCS line is the 
only undetected transition and will be not discussed in the 
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Figure 1. Large-scale view of the surroundings of IRAS 20343+4129 as seen by Spitzer at 8.0 /im. The Cygnus OB2 association is 
located to the south-west as shown by the arrow. The enlargement panel to the right shows the region of our interest as it appears in the 
Spitzer image at 5.8 /im, where we highlight the two infrared sources IRS 1 (saturated in the IRAC image) and IRS 3. In both panels, 
the colour scale units are MJy sr~^. The dashed contours correspond to the 3.6 cm continuum emission detected by Carral et al. 1119991 ) 
with the VLA (D-configuration). 



following. The continuum was derived by averaging the 500 
MHz bands. Visibility data were edited and calibrated with 
the MIRIAD package. A minor flagging of the data was per- 
formed using the UVFLAG task, mainly to remove the inter- 
vals characterized by the bad atmospheric phase coherence. 
The channel spacing and the corresponding 1 a rms noise 
are shown in Table [T] Merging the visibilities obtained in C 
and D conflguration, imaging, deconvolution, and analysis 
of channel maps and continuum were performed using the 
standard tasks of the GILDA^ package (e.g. UVMERGE, 
UVMAP, CLEAN). Images were created applying natural 
weighting to the visibilites. 



2.2 VLA 

120343 was observed with the Very Large Array (VLA|j at 

1.3 and 0.7 cm on 2007 Mar 26 using the array in the D 
configuratiorjj. The phase center of the observations was 
R.A. (J2000) = 20''36'"07'.'51; Dec. (J2000) = 41°40'00''9. 
The data reduction followed the VLA standard guidelines 



1 the GILDAS software is developed at the IRAM 
and the Observatoire de Grenoble, and is available at 
[http://www.iram.fr/IRAMFR/GILDAS 

'' The Very Large Array (VLA) is operated by the National Ra- 
dio Astronomy Observatory (NRAO), a facility of the National 
Science Foundation operated under cooperative agreement by As- 
sociated Universities, Inc. 

^ The VLA continuum observations presented in this work corre- 
spond to project AP525. Another project, AP533, was undergone 
at 3.6, 2 and 0.7 cm, but due to technical problems the data of 
AP533 project were lost. 



for calibration of high frequency data, using the NRAO 
package AIPS. The absolute flux scale was set by observ- 
ing the quasar 1331-1-305 (3C286), for which we adopted a 
flux of 2.52 Jy at 1.3 cm, and 1.45 Jy at 0.7 cm. The quasar 
20154-371, with a bootstrapped flux of 1.39 ± 0.02 Jy at 
1.3 cm and 2.1 ±0.2 Jy at 0.7 cm, was observed regularly to 
calibrate the gains and phases. Final images were produced 
with the robust parameter of Briggs (1995) set to 5, corre- 
sponding to natural weighting. At 7 mm we applied a taper 
at 60 kA with the aim of recovering faint extended emission, 
but no emission was detected at this wavelength. 

The VLA was also used to map the {J,K)={1,1) 
and (2,2) inversion transitions of the ammonia molecule 
on 2001 July 23, with the array in the C conflguration. 
The phase center was set to R.A.(J2000) = 20''36™08':013; 
Dec.(J2000)= -H41°39'56'.'93. The FWHM of the primary 
beam at the observed frequency was ~110", and the range 
of projected baselines was 2.59 to 267.20 kA. The abso- 
lute flux calibration was performed by using 3C286, adopt- 
ing a flux density at 1.3 cm of 2.41 Jy. The phase cali- 
brator was QSO 2013-1-370, with a 1.3 cm bootstrapped 
flux density of 2.34±0.04 Jy, and 3C273 was used as the 
bandpass calibrator. The NH3(1,1) and NH3(2,2) lines were 
observed simultaneously in the 4 IF correlator mode of 
the VLA (with 2 polarizations for each line), providing 63 
channels with a spectral resolution of 0.62 km s~^ across 
a bandwidth of 3.13 MHz, plus an additional continuum 
channel containing the central 75% of the total band- 
width. The bandwidth was centered at the systemic velocity 
VtsR=11.5 km s"^ (Sridharan et al. 2002) for the NH3(1,1) 
line, and at Vlsr=6.5 km s~^ for the NH3(2,2) line (to cover 
the main and one of the satellite components). Data were 
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Table 1. Observed tracers and basic parameters for the CARMA and VLA observations. The synthesised beam and 1 a 
rms for CARMA observations are based on the combined configurations C+D (unless when differently specified) 



Instrument/Tracer ^ 


Frequency _ 


Synth, beam 


Linear 


Av 


Icr rms 


LAS'^ 




(GHz) 


("x") 


resolution (pc) 


(kms-i) 


(Jy beam"-') 


{") 


CARMA 


3 mm continuum 


86.4197 


4.80 X 4.36 


~0.03 


- 


0.0004 


33 


NHaD Jk^.k, = li,i - lo,i (F = 1 - 1) 


85.9263 


3.32 X 2.68 


~0.02 


0.07 


0.045'' 


33 


C2H Nj^p = 13/2,2 — 0i/2,i 


87.3169 


3.35 X 2.72 


~0.02 


0.41 


0.02 


33 


ortho-c-C3ii2 Jk^.k^ = 2i,2 — lo,]^ 


85.3389 


5.81 X 4.54 


~0.034 


0.07 


0.07 


33 


CCS ATj = 76 - 65 


86.1814 


5.67 X 4.84 


~0.035 


0.07 


0.04 


33 


N2H+ Jp^ = I2 - Oi {F = 3^2]f_ 


93.1737 


2.04 X 1.94 


~0.013 


0.065 


0.07 


16 


VLA 


7 mm continuum 


43.3399 


4.67 X 4.57 


~0.03 


- 


0.0004 


18 


1.3 cm continuum 


22.4601 


3.02 X 2.98 


-0.02 


- 


0.000057 


34 


NH3(1,1) 


23.6945 


4.20 X 3.19 


~0.025 


0.62 


0.0035 


34 


NH3 (2,2) 


23.7226 


4.25 X 3.19 


~0.025 


0.62 


0.0035 


34 



"^ For the molecular transitions, in the text we will use the following abbreviations: NH2I3 (li,i — lo,i) = NH2D (1—1); 

C2H (13/2,2-01/2,1) = C2H (1-0); C-C3H2 (2i,2 - lo,i) = C-C3H2 (2-1); N2H+ (I2 ~0i) = N2H+ (1-0); 

" rest frequency of the transition listed in Col. 1; 

^ largest angular scale (at half power) detectable by the interferometer, estimated from the minimum baseline of the array 

configuration, and following the appendix in Palau et al. 1120101) . For the lines observed in C and D configuration, this refers 

to the merged uv coverage; 

sensitivity in the merged C+D channel maps smoothed to a spectral resolution of 0.1 km s"^; 
^ observed in D configuration only; 
^ observed in C configuration only. 



calibrated using standard procedures of AIPS, and imag- 
ing was performed by applying a taper in the -uw-data of 
50 kA and using natural weighting to recover the faint and 
extended emission. The synthesised beams and rms noises 
for all the VLA observations are listed in Table [T] 



3 RESULTS 

3.1 Continuum emission maps 

In Fig.[2]the 3 mm (CARMA, C+D configuration, solid con- 
tours) and 1.3 cm (VLA, D configuration, dashed contours) 
continuum emission maps are shown. As reference, we over- 
plot thel.3 mm continuum emission observed by Palau et 
al. (|2007bl . obtained with about a factor 1.5 better angu- 
lar resolution) as well as the direction of the lobes of the 
^■^CO outflow originated by IRS 1 and the (rough) edge of 
the f an-sha ped H2 emission associated with IRS 3 (Kumar 
et al. l2002D . 

The 3 mm continuum emission is resolved into 5 main 
condensations, which we call MMA, MMB, MMC, MMD 
and MME in order of increasing R.A. The brightest are 
MME and MMA, east and west of IRS 3, respectively. From 
Fig. [2] we note that MME roug hly coi ncides with a 1.3 mm 
continuum peak (Palau et al. l2007a ). while the western- 
most one, MMA, encompasses 3 peaks seen at 1.3 mm with 
SMA. The faintest core, MMC, is detected towards IRS 1. 
Two more 3 mm condensations, MMD and MMB, are not 
detected at 1.3 mm. Specifically, MMB is clearly detected 
close to IRS3, with a shift of only 4" to the northeast. 



The most evident differences among the 3 and 1.3 mm 
conti nuum m aps are that at 3 mm core MME (MM7 in Palau 
et al. l2007bl ) is more extended, and the eastern and west- 
ern cores are connected by a filamentary emission passing 
through IRS 3 which is undetected at 1.3 mm. Both differ- 
ences are probably just the consequence of CARMA being 
sensitive to larger structures than the SMA. We es timated 
that the SMA in the observations of Palau et al. (|2007bl ) 
was sensitive to structures with FWHM < 9" (using the 
mini mum b aseline of the observations and following Palau 
et al. l20ld ). while CARMA using C+D configurations was 
sensitive to structures < 33" (Table [l|, allowing CARMA 
to recover more extended emission. 



As shown in Fig.O the continuum emission at 1.3 cm is 
dominated by one strong and compact source with its emis- 
sion peak (20''36"'07.3', +41°39'52") coincident with IRS 3, 
and matchin g wel l the fan-shaped H2 emission detected by 
Kumar et al. l2002l . A Gaussian fit to this compact centimeter 
source yields a peak intensity of 0.8 mjy/beam, a flux den- 
sity of 1.3 mjy, and a deconvolved size of 4.3" x 1.0", with 
PA=78°, corresponding to a size of 6000 AU in the east- 
west direction (and an unresolved size in the north-south 
direction). In addition to the compact 1.3 cm source asso- 
ciated with IRS 3, there is a secondary peak at around 6a 
which falls 2" to the south of MMD, and faint emission join- 
ing this secondary peak and the centimeter source in IRS 3, 
suggesting that the two peaks of centimeter emission could 
be linked. 
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Figure 2. Continuum emission obtained with CARMA at 3 mm 
(solid contours) and with the VLA at 1.3 cm (dashed contours) 
towards 120343. The solid contours start from the Scr rms level 
(1.2x10^'^ Jy beam~^) and are in steps of Itr rms. The dashed 
first contour and step correspond to the Sc rms level (1.5xl0~^ 
Jy beam~^) and are in steps of la rms The grey scale i ndicates 
the 1.3 mm continuum emission observed by Palau et al. l|2007b|) 
with the SMA (levels = 10% of the maximum, first level = 20% 
of the maximum). The cro sses indicate the two infrared sources 
idenfitifed by Kumar et al. IJ2002V Red and blue arrows highlight 
the direction of the blue - and red-shifted emission detected in 
^^CO (2-1) by Palau et al. l2007bl . The green arc roughly shows the 
fan-shape d feat ure detected in H2 /xm associated with IRS 3 (Ku- 
m ar et al . I2OO2I ). interpreted as an expanding cavity by Palau et 
al. l2007bl . The ellipse at bottom-left corresponds to the CARMA 
synthesised beam (4.8" X 4.4", at P. A. 92°). The dashed ellipse 
at bottom-right is the VLA synthesised beam (3.71" X 3.50", at 
P.A. -56.2°). 



3.2 Distribution of the Integrated intensity of 
molecular line emission 

3.2.1 Molecular tracers observed with CARMA 

The maps of the integrated intensity of the lines observed 
and detected with CARMA (see Table[T} are shown in Fig.[3l 
The emission map of each tracer has been superimposed on 
the images obtained from the Spitzer Space Telescope in the 
four mid-IR IRAC bands (centred at 3.6, 4.5, 5.8 and 8 /im, 
respectively). The location of the near-infrared sources IRS 
1 and IRS 3 are also indicated, as well as the direction of the 
outflow centred on IRS 1 and the edge of the H2 emission 
associated with IRS 3, as in Fig. [2] We also superimpose 
the 1.3 cm continuum emission detected in this work (see 
Fig- 01 which marks clearly the Hll region associated with 
IRS 3. 

The molecular gas seems to be squeezed in between 
the two dominant mid-IR sources IRS 1 and IRS 3 in all 
molecules except N2H"'". Other compact mid-IR sources are 
present in the region but do not seem to be associated with 
any clear molecular counterpart. The diffuse IR nebulosity, 
especially evident in the 5.8 and 8 ^m bands, is probably 
emission from small dust grains distributed around 120343 
becoming brighter at longer wavelengths. Some of the dif- 



fuse emission dete cted a t 8 fj,in may also be PAHs emission 
(e.g. Peeters et al. l2002l ). 

The morphology of the integrated intensity of C-C3H2 
(2-1) and C2H (1-0) delineates clearly a cavity around IRS 3 
(top panels of Fig. [Sjl , providing a strong support to the hy- 
pothesis proposed by Palau et al. l|2007bl ). namely that IRS3 
is most likely a more evolved intermediate-mass star creating 
a cavity. The C2H emission is more extended than that of c- 
C3H2, perhaps due to the smaller sensitivity that we have in 
the C-C3H2 channel maps (see Table [TJ. Specifically, a nar- 
row filament extended ~ 30" is clearly detected north-east 
of the field center (Fig. |31 top-left panel), inclined roughly 
as the tail of the mid-IR cometary shape (see Fig. 1), sug- 
gesting that the two features can have the same origin. 

The bottom panels in Fig. |3] show the integrated emis- 
sion of the two Nitrogen-bearing species N2H+ and NH2D. 
The emission in N2H"'" consists mainly of one cloud to the 
east of IRS 3 elongated in the southeast-northwest direc- 
tion, and extending up to IRS 1, and two smaller clouds, 
one immediately to the south-west of IRS 3 with no con- 
tinuum counterpart (called IRS3-SW), and the other asso- 
ciated with MMA. In addition, there is one clump about 
1' to the west of IRS 3, almost at the border of the bright 
rim, which falls on a region with no infrared emission asso- 
ciated. On the other hand, the emission of NH2D consists 
mainly of one filament elongated in the east- west direction, 
passing through IRS 3, and with some emission at IRS3- 
SW. The N2H+ clump 1' to the west of IRS 3 is detected 
also in NH2D, but looks more extended in NH2D. However, 
this can be just an efl'ect of the different angular resolution 
and different filtering of extended emission, as NH2D was 
observed with C+T) configuration while N2H"'' was observed 
in C configuration only (see Table [T}. 

If we put together the two mostly extended molecu- 
lar tracers, namely C2H and NH2D, we can notice a sort 
of 'snake-like' filament of molecular gas (Fig. 3} extending 
from the south-western side of 120343, clearly detected in 
NH2D, up to the north-eastern corner, in which a long and 
narrow filament is detected in C2H. The bulk of the emis- 
sion is in between IRS 1 and IRS 3. This 'snake-like' filament 
matches very well the 1.2 mm continuum emission detected 
with MAMBO (Beuther et al. 2002, Fig. ID), and its SW- 
NE inclination follows roughly the 'head-tail' orientation of 
the mid-IR diffuse cometary emission (Fig. [T|, suggesting a 
possible common origin. However, the highest sensitivity re- 
gion of the CARMA maps, i.e. the field of view of the 10 m 
dishes (~ 73"), includes only the central region of the fil- 
ament (see Fig. 3}. For this reason, in this work we focus 
on the centre of 120343, where the interaction between the 
two brightest IR sources and the surrounding molecular gas 
seems predominant, and plan a large interferometric mosaic 
which will allow us to unveil the overall distribution of the 
molecular gas and its relation with all the IR sources. 



3.2.2 Ammonia (1,1) and (2,2) inversion transitions 

The integrated intensity maps of NH3 (1,1) and (2,2) are pre- 
sented in Fig. [5] In the (1,1) line, the emission resembles that 
seen with CARMA in the 3 mm continuum. Four main peaks 
are detected, which roughly correspond to those detected 
in the 3 mm continuum. On the other hand, none of the 
1.3 mm continuum peaks identified by Palau et al. l|2007bl ) 
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Figure 3. Integrated maps of the molecular lines detected with CARMA (Table [T]| towards 120343 and superimposed on the images 
obtained in the Spitzer-IRAC bands (3.6 /xm, 4.5 /im, 5.8 lira, 8 /im, colour scale in MJy sr~^ units). In each panel, black contours 
depict the velocity-averaged emission of the transition labelled at the top-right corner (the corresponding synthesised beam is shown 
in the bottom-right corner). For all lines, the emission has been averaged over all the velocity channels with signal, except for N2H+ 
(bottom-left panel) for which the integrated emission was averaged over the main group of hyperfine components only. For C2H and 
NH2D, first contour and step correspond to the 20% level of the maximum (corresponding roughly to the 3(7 rms level of the averaged 
map), while for C-C3H2 and N2H+ contours start from the 30% level of the maximum and are in steps of 20%. The white contours 
represent the VLA 1.3 cm (K-band) continuum image (same contours as in Fig. [2]|. The position of IRS 1 and IRS 3, the associated 
outflow and fan-shaped H2 emission are shown as in Fig. [2] 



exactly coincides with any of the NH3 (1,1) peaks, despite 
the similar angular resolution. The NH3 (1,1) emission re- 
veals two main clouds, one to the east and the other to the 
west of IRS 3. The eastern cloud includes MMD and MME 
and appears elongated in the southeast-northwest direction, 
similar to the N2H^ eastern cloud and to the 3 mm contin- 
uum emission. The western cloud has a condensation near 
MMB (near IRS 3) and another condensation near MMA. 
The overall emission in the western cloud is elongated in 
the east-west direction. We stress that NII3 (1,1) emission is 
marginally detected also towards IRS 3 and the cavity driven 
by it, as for the 3 mm continuum. On the other hand, the 
millimeter continuum sources associated with IRS 1 as well 
as the ^■^CO outflow lobes are not detected in ammonia. The 
NII3 (2,2) transition is clearly detected towards the eastern 
side of MME, where 2 peaks are resolved. The NH3(2,2) 
line emission around IRS 3 and towards the west is clumpy, 
peaking towards MMA and near MMB and IRS3-SW. 



From a comparison of the N-bearing to the C-bearing 
species, both C2H and C-C3H2 highlight clearly the cavity 
associated with IRS 3 (see top panels in Fig.|3J, while NH2D 
and NH3 trace emission extending east- west passing through 
IRS 3 (see Fig. [5] and bottom panels of Fig. |3} more similar 
to that seen in the 3 mm continuum map. Thus, there seems 
to be a chemical dichotomy in 120343 among Carbon- and 
Nitrogen-bearing molecules. The exception is represented by 
N2H"'", detected away from IRS 3. This molecule appears to 
trace the part of the cloud less disrupted by the expanding 
cavity. 

The spectra of the NH3(1,1) and NH3(2,2) emission, in- 
tegrated within the 5cr contour polygon of the 3 mm contin- 
uum image (except for MMB and MMD, where we used the 
3(7 contour) are shown in Fig. (6] For IRS3-SW, we used the 
5(7 rms contour of the NH2D (1-1) line integrated emission 
(Fig- El bottom-right panel). To compare the different molec- 
ular species, in Fig.|6]we also show the integrated spectra of 
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Figure 4. Integrated emission of molecular tracers (white con- 
tours) superimposed on the image of 120343 in MAMBO at 1.2 
mm (Beuther et al. 2002). The solid white contour depicts the 
15% level of the zero-order moment map of C2H (1—0). The 
dashed white contour corresponds to the 15% level of the zero- 
order moment map for the NH2D (1—1) line. The ellipse in the 
bottom-left corner shows the CARMA synthesised beam for the 
C2H map and it is almost coincident to that of the NH2D map. 
The grey dotted-dashed circles represent the primary beams of 
the CARMA 10 and 6 m dishes (73 and 121 "), respectively. 



NH2D and N2H extracted using the same polygons. We do 
not show the spectra towards MMC because this core is un- 
detected in NH2D and NH3 (2,2), and marginally detected 
in the other lines. Among all the NH3 (1,1) spectra, MMA 
shows the strongest emission, and MME shows the broad- 
est lines, of up to 2.2 km s~ . Such a large line broadening 
in MME could be due to a double velocity component, as 
suggested by the N2H"'" spectrum which has 10 times better 
spectral resolution. Concerning MMB, MMD, and IRS3-SW 
NH3 (1,1) spectra, it is striking the anomaly seen in the in- 
ner satellite hyperfine lines, with one inner satellite clearly 
detected above 5cr and the other inner satellite remaining 
undetected. The anomaly for the non-LTE case due to hy- 
perfine selective photon trapping affects only the oute r sate l- 
lites (red stronger than blue, Stutzki & Winnewisser [l985|) , 
allowing us to rule out this possibility in 120343. Rather, 
anomalies of one inner satellite being stronger than the other 
have been observ ed in several works (Lee et al. 2002, Long- 
more et al. |2007| . Purcell et al. |2009| ) and explained as being 
due to s ystem atic mo tions, following the theoretical work 
of Park (|200ll ). Park (|200ll ) shows that if the core is con- 
tracting the inner blue satellite should be stronger than the 
inner red satellite for a systematic motion in the range of 
0.4-1 km s~^, and for a range of H2 number densities and 
NH3 column densities which are consistent with those de- 
rived by us (as we will show in Sect. 13.4]) . This is the case of 
MMD and IRS3-SW. On the other hand, if the core is under- 
going expanding motions, the prediction is that the inner red 
satellite should be stronger than the inner blue satellite, as 
seen for the case of MMB. Thus, it seems that for these three 
clumps the NH3(1,1) anomalous intensity of the hyperfine 



components is consistent with contracting/expanding mo- 
tions. 

The MMA NH3(1,1) spectrum is a very special case, 
as the inner satellites are detected at an intensity smaller 
than that expected for LTE conditions (maximum main-to- 
satellite ratio in LTE is 3.6, while the ratio for MMA is 4.3). 
This could be explained if the opacity is high and the ex- 
citation temperature for the main line and the satellites is 
different, with the satellites having a lower excitation tem- 
perature (and hence a smaller main beam temperature). A 
detailed discussion on the temperature ratio between main 
line and satellites leads to a possible non isothermal core 
made of two layers, with the external one being hotter than 
the inner one. If this interpretation is correct, core MMA 
could be heated externally, perhaps by IRS 3 and/or the in- 
frared sources west of 120343 (see right panel in Fig. [T] and 
Fig- El). We give details of this explanation in Appendix A. 

We fitted the NH3, NII2D, and N2H+ spectra in order 
to derive the physical parameters of the gas traced by these 
molecules. To take into account the line hyperfine structure, 
we followed the method described in the CLASS user man- 
uaB The NH3 (1,1), NH2D (1-1) and N2H+ (1-0) lines were 
fitted this way, while we fitted the NII3 (2,2) lines with Gaus- 
sians. The derived fit parameters are reported in Table (2] 
except the line velocities that will be extensively discussed 
in Sect. 13.31 These parameters have been used to derive the 
molecular column densities, the derivation of which will be 
described in Sect. 14.21 To fit the hyperfine structure of MMB, 
MMD, and IRS3-SW, we used only the detected satellite and 
the main line, in order to obtain reliable opacities. For NH3, 
in Table [5] we list also the rotation temperature, Trot , ob- 
tained from the (2,2)-to-(l,l) intens ity ra tio following the 
method outlined in Busquet et al. (|2009l ) . which is based 
on the discussion presented in Ho & Townes (1983). These 
range from 13 K (in MMD) to 23 K (in IRS3-SW). 

In order to obtain a first approach to temperature vari- 
ations of the dense gas across 120343, we computed the ratio 
of the integrated NH3 (2,2)/(l,l), which has been shown to 
be a reasonable approach to the kinetic temperature (e.g.. 
Torrefies et al. 1993; Zhang et al. 2002). The result is shown 
in Fig. [7] As can be seen from the figure, the ratio is largest 
towards three main regions: near IRS 3 (MMB, and IRS 3- 
SW), to the north of MMA, and to the eastern side of MME. 
Thus, the spectra toward MMA, which can be explained 
with external heating, is consistent with the 22/11 ratio map 
which reveals that the heating comes most likely from the 
north of MMA. Interestingly, the north of MMA is spatially 
coinciding with the redshifted CO outfiow lobe driven by 
IRS 1, suggesting that the outflow and the dense gas are 
interacting (Fig [7}. Most intriguing is the high ratio seen at 
the eastern side of MME, which extends all along the north- 
south direction and is seen westwards of the H2 extended 
emission reported by Kumar et al. (2002). Such a spatial 
coincidence is suggestive of a possible relation between the 
high 22/11 ratio and the H2 emission. Finally, the high 22/11 
ratio near IRS 3 could be indicative of direct heating by the 
UC Hll region associated with IRS 3. 



* The CLASS program is part of the GILDAS software, devel- 
oped at the IRAM and the Observatoire de Grenoble, and is avail- 
able at htt p://www.iram.fr/IRAMFR/GILDASI 
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Table 2. NH3, N2H"'", and NH2D line parameters for the 3 mm continuum cores (except MMC) and the molecular core 
IRS3-SW. The parameters have been derived from the spectra shown in Fig. |6] using the fitting procedure described in 
Sect. 13X21 
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23 


0.72 
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" = Total opacity of the line; 

' = in this core the excitation temperature for the main and the satellite lines of the NH3 (1,1) transition is probably 

different. See Appendix A for the detailed discussion on the derivation of the physical parameters from NH3. 
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Figure 6. Spectra of the detected transitions of NH3, NH2D and N2H"'" at the position of the 3 mm continuum condensations (except 
MMC, undetected in all the lines) and the molecular condensation IRS3-SW. The polygon used to extract the spectra corresponds to: 
the 5cr rms contour of the 3 mm continuum image for MMA and MME; the 3a" rms contour of the 3 mm continuum image for MMB and 
MMD; the Scr rms contour of the integrated emission of NH2D (1-1) for IRS3-SW. 
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Figure 5. Top panel: integrated emission map of the ammonia 
(1,1) inversion transition (contours) observed with the VLA to- 
wards 120343, superimposed on the 3 mm continuum (CARMA, 
grey-scale). First contour level and step are the 10% of the peak 
(224 Jy beam~^ m s~^). The ellipse at bottom-left represents 
the VLA synthesised beam (4f'20 X 3f' 19, at P.A. ~ 11°). Purple 
triangles pinpoint the 3 mm continuum peaks (see Fig. [2| . All 
other symbols (crosses, arrows and semiellipse) are the same as 
in Fig. [2] 

Bottom panel: same as top panel for the (2,2) inversion tran- 
sition. First contour level and step are the 15% of the peak 
(57 Jy beam~^ m s~^). The synthesised beam in this case is 
4f'25 X 3f'19, at P.A. ~ 13°. 



3.3 Kinematics 

To inspect the velocity field, we have extracted from the 
interferometric channel maps spectra of the molecular tran- 
sitions detected on grids with regular spacings (1.5" x 1.5" 
for C2H, 2.5" X 2.5" for C-C3H2, 1" x 1" for N2H+, 1.6" x 
1.5" for NH2D and NH3). The spectra extracted have been 
fitted following different methods: for C2H (1-0), C-C3H2 
(2-1) and NH3 (2,2) we assumed Gaussian lines, while for 
N2H+ (1-0), NH2D (1-1) and NH3 (1,1) we fitted the lines 
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Figure 7. Colorscale: ratio of the integrated NH3 emission 
(2,2)/(l,l). Red triangles indicate the 3 mm continuum peaks, 
which correspond to (from east to west) MME, MMD, MMC, 
MMB and MMA. The black-dashed contours indicate the VLA 
1.3 cm emission (2.5, 5, 7.5, 12, 15 times 50 /ijy beam— ^). 
The blue/red solid contours depict the high-velocity CO (2—1) 
blueshifted/redshifted emission (Palau et al. |2007b|). Blue con- 
tours range from 10 to 99% in steps of 12% of the peak intensity 
(21.335 Jy beam"'^ km s"'^); red contours are 6, 12 to 99% in 
steps of 12% of the peak intensity (41.781 Jy beam~^ km s~^). 



following the method described in Sect. 13.2.21 to take into 
account their hyperfine structure. 

3.3.1 Line peak velocities 

Fig. [8] shows maps obtained from the line peak velocities. 
In all tracers the radial velocities are predominantly blue- 
shifted to the west and red-shifted to the north-east. The 
east-west velocity gradient is not uniform and suggests a 
possible torsion of the gas. Interestingly, the inclination of 
this gradient with respect to the line of sight is opposite to 
that of the outfiow associated with IRS 1, since the blue lobe 
of the outflow is located on the side where the dense gas is 
red-shifted, and vice- versa. In general, all tracers with emis- 
sion near IRS 3 (in MMB) show that the gas is blueshifted 
at this position. 

3.3.2 Line widths 

Maps of the line widths are presented in Fig. [9] The mea- 
sured line broadenings are generally a factor > 3 higher 
than the thermal broadening, expected to be of the or- 
der of ~ 0.1 — 0.3 km s~^, indicating that the gas kine- 
matics is largely dominated by non-thermal motions. This 
finding confirms previous studies in simil ar inter mediate- to 
hi gh-ma ss protocluste rs (e.g . Palau et al. l2007al . Fontani et 
al. l2009l . Busquet et al. l20ld ). and represents one of the most 
important differences between these dense cores and those 
observed in low-mass star forming regions, where the line 
widths are dominated by thermal broadening, even in clus- 
tered environ ments (e.g. Kirk et al. |2007| . Walsh et al. |2007| , 
Bourke et al. l201lh . 

In summary we highlight three regions where the gas 
is more turbulent (Fig. (Qjl: (i) around IRS 3, especially in 
between IRS 1 and IRS 3 (see the C2H and C-C3II2 line 
widths in Fig. [5|; (ii) north of MMA, near the red lobe; 
(iii) towards MMA. The turbulence enhancement around 
IRS 3 is easily explained by the expading cavity, while in 
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MMA and north of it could rather be due to the red-lobe 
of the outflow. On the other hand, we stress that towards 
MME and in between IRS 1 and IRS 3 some spectra show 
a double- velocity component (see e.g. the isolated hyperfine 
component in the N2H+ (1-0) line towards MME in Fig. [6]). 
These are especially evident in the N2H"'', NH3 and NH2D 
spectra, but the second component can be seen also in some 
spectra of C2H and C-C3H2 close to the cavity. Therefore the 
large broadening above IRS 3 and in MME could just be due 
to the superposition of two unresolved velocity components. 
Interestingly, for C2H and N2II"'" the line broadening is 
relatively small close to the blue lobe of the outflow associ- 
ated with IRS 1 and larger besides the red lobe, while the op- 
posite seems to occur for the C-C3II2 (2-1) line (top-middle 
panel in Fig. [9l). This could indicate a selective influence of 
the flow on the different molecules in the surrounding ma- 
terial. However, we stress that the C-C3II2 emission is very 
faint at the borders of the region plotted in Fig. [§1 where 
the fit results are affected by large relative errors. 



3.3.3 Position-velocity plots 

In order to further study the velocity field of the NII3 emis- 
sion near IRS 3, we computed position- velocity plots for the 
NII3 (1,1) line in the east-west direction, and centred at off- 
set 0",— 2" with respect to the CARMA phase centre. The 
final plot is shown in panel "a" of Fig. [TO] The NH3 (1,1) 
emission in the position-velocity plot shows two main peaks, 
one corresponding to the eastern cloud and the other cor- 
responding to the western cloud, and both peaks are linked 
through fainter emission which overall shows a U-like struc- 
ture. Such a feature resem bles the shape predicted by the 
model of Arce et al. (|201in for an expanding bubble. In this 
model, an expanding shell would appear as a ring in the (p- 
v) plots (see their Fig. 5), while we only see the blue-shifted 
half of it. However, a U-like feature can be explained if the 
source driving the bubble is slightly displaced (behind the 
bubble) with respect to the surrounding molecular gas, so 
that we mainly see the gas which is moving towards us. In 
fact, the tails of the U-like feature are found, as expected, 
at approximately the systemic velocity. This suggests that 
IRS 3 may be pushing the surrounding dense material away 
(either through its winds/radiation or through the associ- 
ated Hll region), with an expansion velocity (difference be- 
tween the 'tip' of the U-like feature and the 'tails') of about 2 
km s^^ (see panel "a" in Fig. llOp . A similar expanding shell 
was recently found around an infrared-source at the centre of 
a region devoid of gas emission in the intermediate- to hig h- 
mass protocluster IRAS 053454-3157 (Fontani et al. 120121 '). 



3.4 Physical parameters from 3 mm and 1.3 cm 
continuum emission 

3.4.1 3 mm 

We have identified 5 main condensations in the 3 mm con- 
tinuum emission map (see Sect. [5TT1) . Their peak positions 
are given in Cols. 2 and 3 of Table[3l We also list the angular 
{9) and linear (D) diameters (Cols. 6 and 7 in Table[3]) com- 
puted assuming the sources are Gaussians, and deconvolving 
the contour at half maximum with a Gaussian beam with 
HPBW corresponding to the geometric mean of the minor 
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Figure 10. a) NH3 (1,1) position-velocity plot centred at offset 
(0,-2") with respect to the CARMA phase centre and along a cut 
in the east-west direction. The NH3 (1,1) data cube was convolved 
with a beam of 5"x2" and PA=0° to increase signal-to- noise. 
Contours are: —3, 3, 6, 9, 15, 21, and 27 times 0.001 Jy beam~^. 
b) NH3 (1,1) integrated emission map (same contours as in Fig. 5) 
with the cut of the position- velocity plot of panel "a" indicated 
by a dashed line. Contours start at 12% of the peak intensity, 
224 Jy beam~^ m s~^, and increase in steps of 10%. The VLA 
synthesised beam is shown in the bottom-right corner. 



and major axes of the CARMA synthesised beam (see Ta- 
ble[T|. Because the contours at half maximum were blended, 
at the edge between two cores we decided to separate the 
emission arising from the different condensations identify- 
ing the peaks and considering the first unblended contour. 
The same criterion was applied to derive the integrated flux 
density, Fv, given in Col. 5 of Tabled 

From Fv, we have computed the mass of the conden- 
sations under the assumptions that the dust millimeter- 
continuum emission is optically thin, and that the dust tem- 
perature equals the gas kinetic temperature. This latter hy- 
pothesis implies coupling between gas and dust, which is 
a reasonable assumption for H2 volume densities above 10^ 
cm~"^. Under these assumptions, the gas mass can be derived 
from the formula: 



M = 



F.d^ 



B,{T)K 



(1) 



where d is the source distance, B^{T) is the Planck func- 
tion at dust temperature T, and fc^ is the dust opacity per 
unit dust mass. For this latter, we e xtrapo lated the value 
at 230 GHz given by Kramer et al. (|l998l '). ^230 = 0.005 
cm^ g~^ (which assumes a gas-to-dust ratio of 100), through 
the power-law k^, = A;23o[i'(GHz)/230]'^. We have assumed 
13 — 2, which is a typical value derived for dusty envelopes of 
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Figure 8. Maps of the line peak velocity for all detected molecular lines. Top panels show (from left to right): C2H (1—0), C-C3H2 (2—1) 
and N2H+ (1-0). Contours are in steps of 1 km s~^, and range from: 8 to 11 km s~^ for C2H; 10.5 to 13.5 for C-C3H2; 9 to 12 km s^^ 
for N2H+. Bottom panels show (from left to right): NH3 (1,1) NH3 (2,2) and NH2D (1-1). Contours are in steps of 1 km s~^, and 
range from: 9.5 to 13.5 km s"! for NH3 (1,1); 9 to 13 km s"! for NH3 (2,2); 8.5 to 12.5 km s'^ for NH2D. In each panel, the ellipse 
in the bottom-right corner represents the CARMA or VLA synthesised beam. The black triangles pinpoint the 3 mm continuum peaks 
(Sect. [3T1 . The position of IRS 1 and IRS 3, as well as the associated outflow and fan-shaped H2 emission, are depicted in each panel 
as in Fig. [2] 
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Figure 9. Same as Fig. [8]for the lino widths. Contours are in steps of 0.5 km s ^ and range from: 0.5 to 3 km s ^ for C2H 
km s-i for C-C3H2; 0.4 to 2.4 km s"! for N2H+; 0.6 to 2.6 km s"! for NH3 (1,1); 0.8 to 3.2 km s'^ for NH3 (2,2) ; 0.3 to 1 
for NH2D (1-1). 
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emission measure: 3.4 10 cm pc 
electron density: 1.0 10 cm 
' ionizing photons: 2.0 lO'^'^s"^ (B2 ZAMS) 

spectral index: 0.1 +/-0.2 
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Figure 11. Spectral energy distribution of IRS 3 in the cm/mm 
range. The physical parameters have been derived assuming the 
emission at 1.3 cm comes from an optically thin Hll region. 



massi ve (p roto) stellar objects (e.g. Molinari et al. 2000, Hill 
et al. 120061 ). As gas temperature, we have taken the kinetic 
temperature obtained by extrapolating the rotation temper- 
ature derived from the ammonia (2,2)/(l, 1) line ratio for 
each core (see Sect. 13. 2T2] and Table[2]) following the empiri- 
cal approximation method described in Tafalla et al. (|2004r ). 
The kinetic temperatures derived this way are listed in Col. 4 
of Table m and are in between 14 and 25 K, higher (a factor 
1.5 - 2) than the values measured ty picall y in low-mas s clus - 
tered starless cores (e.g. Andre et al. 120071 . Foster et al.|2009|) 
which are around 10-13 K. For MMC, undetected in NH3 
(2,2), we decided to give a range of masses computed in the 
temperature interval 15-30 K. 

The resulting masses are listed in Col. 8 of Table |31 
All fragments have masses consistent with intermediate-to 
high-mass embedded objects. The most massive one is MME 
(23 Mq). By assuming spherical and homogeneous cores, 
we have derived the average H2 volume and column den- 
sities. The average volume and column densities (given in 
Cols. 9 and 10 of Tabled are of the order of lO'^-'^cm"^ and 
10^^~^'*cm~^. Such high column densities are consistent with 
being the birthplaces of interm ediate - and/or high-mass ob- 
jects (e.g. Krumholz & McKee bOOi ). 



34.2 1.3 cm 



In Sect. l3.l1 we showed that the centimeter emission is domi- 
nated by one single source associated with IRS 3 and extend- 
ing towards the east. In addition, a secondary source near 
MMD was also identified. A 3.6 and 6 cm source assoc iated 
with IRS 3 is already reported by Miralles et al. (jlQOJ) and 
Carral et al. (1999). However, from these two measurements 
only, and taking into account the uncertainties, the spectral 
index of the centimeter source associated with IRS 3 could 
not be well determined (e.g., Palau et al. l2007bi l . Our new 
measurement at 1.3 cm allows to better constrain the spec- 
tral index of the source to 0.1±0.2 (Fig. 1111) . Such a spectral 



index is consistent with optically thin free-free emission fa- 
voring the interpretation that the centimetre emission comes 
from an Hll region rather than a thermal radio jet (which 
typically have steeper spectral indices). We calculated the 
physical parameters of the ionised region at 1.3 cm assum- 
ing the emission is optically thin, and obtained an emission 
measure of 3.4 x lO"* cm~^ pc, characteristic of UC Hll re- 
gions, and a flux of ionising photons of 2 x 10** s~^, consis- 
tent with the Hll region being ionised by an early-type B2 
star (Panagia 1973). Interestingly, the extension of such an 
UC Hll region towards the east is similar (although larger in 
size) t o the extension seen in [Nell] emission by Campbell 
et al. (I2OO3), and which is interpreted as due to the expan- 
sion of the ionised gas and disruption of the natal envelope. 
In this context, the secondary centimetre peak near MMD 
could be related as well to the expanding ionised gas. We 
estimated the possible contribution of free-free emission to 
the flux measured at 3 mm for MMB, and is of ~ 0.3 mjy 
(4cr of 1.3 cm observations, and using the spectral index of 
0.1), out of 4.6 mJy, or 6%. Thus, thermal dust emission is 
the main contribution to the 3 mm continuum emission in 
MMB. Finally, we estimated an upper limit for the 3 mm 
emission associated with IRS 3 of 2.3 mJy, measured as the 
3 mm flux density inside the 4(t contour of the centimeter 
emission (Fig. llip . and we cannot rule out the possibility of 
the UC Hll region being still associated with remnant natal 
dust, although a projection effect could be also possible. 



4 DISCUSSION 

4.1 Column densities of the PDR tracers C2H 
and C-C3H2 

The two carbon-bearing species C2H and C-C3H2 are among 
the most abundant simple carbon-chain molecules detected 
in the interstellar medium, and ar e believed t o be g ood trac- 
ers ofPDRs (Lucas & Liszt I2OOOI . Pety et al. |2005| . Gerin et 
al. |201l|). C2H is formed either from photodissociation of 
acetylene (C2H2) followe d by d issociative recombination of 
C2HJ (Mul & McGowan [19801 ) or through neutr al-neu tral 
reaction between C and CH2 in hot gas (Sakai et al.|2010|). c- 
C3H2 is believed to be formed by dissociative recombination 
of c-CsW^ . Both species benefit from the presence of atomic 
carbon not locked in CO, and a good correlation between 
the two tracers has been found at the ill uminated sur face o f 
the Horsehead nebula (Pety et al. |2005| . Gerin et al. boogj ). 
as w ell as in both diffu se an d translucent clouds (Lucas & 
Liszt 2000, Gerin et al. l20lj ). 

We have investigated the relation among the two species 
in 120343. For this purpose, we have extracted the spectra 
of C2H (1-0) and C-C3H2 (2-1) on a grid of spacing 2.5" 
X 2.5" (roughly half of the CARMA synthesised beam at 
the frequency of the C-C3H2 (2-1) transition), and fitted 
the spectra with Gaussian lines. Then, from the integrated 
intensity obtained from the fits, we have computed the col- 
umn densities assuming that both lines are optically thin. 
This assumption is mandatory because the opacity of the 
lines cannot be directly measured (for the C2H (1-0) line 
we observed only the main hyperfine component, and we do 
not have isotopologues for C-C3H2). However, given that the 
line profiles generally do not show effects due to high optical 
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Table 3. Peak position, angular and linear diameter, integrated flux density, mass, H2 volume and column density of the 
3 millimeter condensations detected with CARMA. The masses are computed for /3 = 2, and assuming the temperatures 
derived from ammonia. The H2 volume and column densities are calculated assuming a spherical source with diameter equal 
to the deconvolved level at half maximum. F^ has been obtained by integrating the continuum flux density inside the So- 
contour level. 



Core 


Peak 


position 


Tu " 


F, 


Diameter 


-^^cont 


"H2 


N(B.2) 




R.A.(J2000) 


Dec.(J2000) 






^s 


D 










20** 36™ 




(K) 


(mjy) 


(") 


(pc) 


(M0) 


(xlO'5 cm-3) 


{xl023 cm-2) 


MMA 


06^37 


+41°39'59f'l 


22'' 


7.4 


7.3 


0.049 


21 


6.7 


10.3 


MMB 


07^54 


+41°39'55"6 


25 


4.6 


6.8 


0.046 


11 


4.5 


6.4 


MMC 


07^72 


+41°40'07"2 


- 


1.9 


5.5 


0.037 


8c^4d 


6.2^2.9'' 


7.1^3.3'^ 


MMD 


08^23 


+41°40'02f'0 


14 


3.7 


6.3 


0.043 


17 


8.7 


11.4 


MME 


08^26 


+41°39'54f'7 


22 


8.3 


8.0 


0.054 


23 


5.7 


9.6 



" = Derived from Trot as explained in Sect. 13.41 

^ = see Appendix A for the derivation of Tj^ from ammonia for this core; 

" = assuming a temperature of 15 K; 

'^ = assuming a temperature of 30 K. 



depths, we are confident that the assumption is reasonable. 
We used the general formula for optical ly th in transitions 
(compare to, e.g., Eq. A. 3 in Pillai et al. l2007h : 



iVtot = 



3h Q{T,^) 



W 



(^) 



87r3 Sfi^ J^Tex)- J.(rBG)eWfcT,. 



(2) 



where: Ej and S are energy of the upper level and line 
strength, respectively, W is the integrated intensity of the 
line, Q(rex) is the partition function at the temperature 
Tex, I' the line rest frequency, Jv{Tcx) and J,^{Tbg) are 
the equivalent Rayleigh- Jeans temperature at frequency f 
computed for the excitation and background temperature 
(Tbg ~ 2.7 K), respectively; n the molecule's dipole mo- 
ment (0.77 Debye for C2H and 3.27 Debye for C-C3H2). For 
C2H, in Eq. ^ W has been obtained by multiplying the 
integrated emission of the hyerfine component observed for 
its relative intensity (0.416). As excitation temperature, we 
have assumed a reasonable value of 20 K based on the ex- 
citation temperatures computed for the other lines (see Ta- 
ble [2]). The values of E,j, Sfj,^ and Q have been taken from 
the Cologne Database for Molecular Spectroscopy (CDMqj, 
Miiller et al. 120011 '). For this latter, we have extrapolated the 
values tabulated to an excitation temperature of 20 K. For 
C-C3H2, the ortho-/para- ratio is included in the partition 
function. 

The results are shown in Fig. 1121 The column density of 
C2H is of the order of 10^* cm~^ across the cloud (Fig. 1121 
top panel), while that of C-C3H2 is of the order of 10^^ cm^^ 
(Fig. ll2l middle panel). The C2H column densities are gener- 
ally larger than those found by Gerin et al. (i2011i ), who mea- 
sured column densities of 10^^ cm~^, while those of C-C3H2 
are more consistent. Their ratio is on average of the order of 
200 - 400 (Fig. 1121 bottom panel), i.e. one order of magni- 
tude larger than th e valu e 20 - 30 measured in translucent 
clouds (Gerin et al. l201l|). as well as in diffuse high latitude 
clouds (Luc as fc Liszt [2000|) and in the Horsehead nebula 
(Pety et al. [2005|). On th e othe r hand, the chemical models 
of PDRs in Gerin et al. (|201ll ) seem to be more consistent 
with our observational results rather than with theirs, be- 
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cause the models predict total column densities consistent 
with our values for both species, and ratios of the order of 
100 or even more (see their Table 5). Interestingly, we find 
the largest ratios (around 400-800) close to the outflow lobes 
and to the east and west of the cavity walls, where the gas 
is probably less dense because most disrupted. Significant 
enhancement can be noticed also in the eastern clump, in 
between MMD and MME, where 1.3 cm emission is detected 
(see Sect. l3.l1 and Fig. [21). This would confirm strongly that 
both molecules are produced in PDR regions, and that they 
are maybe tracing a low density envelope in which the dense 
cores detected in N2H+ and NH3 are embedded. 



4.2 Chemical differentiation and nature of the 
3 mm continuum cores 

Studies of intermediate- and high-mass star-forming regions 
suggest that the relative abundance ratio NH3-to-N2H"'" is 
an evolu tionary indic ator f or dense cores (e.g. Palau et 
al. l2007al . Fontani et al. l2012h : cores with no signs of star for- 
mation typically have larger NH3-to-N2H+ column density 
ratio than cores associated with active sta r form ation. Also, 
Fontani et al. (2008) and Busquet et al. l|20ld ) have mea- 
sured that the deuterated fraction (i.e. the column density 
ratio of a deuterated specied to that of the hydrogenated 
counterpart) of NH3 and N2H^ is of the order of 0.1 in 
pre-protostellar core candidates, as high as in low-mass pre- 
stellar cores, while it is lower in more evolved objects. With 
this in mind, we have investigated the evolutionary stage 
of the millimeter cores in 120343 based on the column den- 
sity ratios NH3-to-N2H+ and NH2D-to-NH3. From the line 
parameters derived in Sect. l3.2T2l fsee Table[2]) we have com- 
puted the total column densitie s of N H3, NH2D and N2H'*' 
from Eq. (Al) of Caselli et al. (|2002l ') because all lines are 
optically thick (except NH2D (1-1) in MMD, but the opac- 
ity is well-constrained). As for the other parameters, for a 
detailed discussion on the NH3 data of MMA see the Ap- 
pendix A. 

In Table[4]we report the column densities of NH3, NH2D 
and N2H"'", and the column density ratios NH2D-to-NIl3 
and NH3-to-N2H+. The NH2D-to-NH3 ratio is of order 0.1 
(from ~ 0.07 to 0.15), and does not change greatly from 
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Table 4. Total column densities of NH3, NH2D and N2H"'" for the 3 mm continuum cores, except MMC, and the additional core 
IRS3-SW. 



A^(NH3) 



iV(N2H+) 



Af(NH2D) JV(NH2D)-to-Af(NH3) Ar(NH3)-to-Af(N2H+) 



(xlO^s cm-2) (xlQi'^ cm-2) {xW^ crci~^) 



MMA 


> 2.41'^ 


MMB 


2.3(0.6) 


MMD 


3.1(0.7) 


MME 


4.4(0.7) 


IRS3-SW 


1.8(0.4) 



7(2) 

13(4) 

5.5(0.7) 

6(3) 



2.2(0.2) 
2.6(0.3) 
2.1(0.2) 
6.8(0.7) 
2.1(0.2) 



< 0.09 

0.11(0.04) 

0.068(0.002) 

0.15(0.04) 

0.12(0.04) 



>34 

24(13) 
80(23) 
30(22) 



N(C2H)' -I J cm 








20 







A R.A 


(arcsec) 


-20 


10 




M 


P 


l« 




' NCc-CjHj) y 





- 


1 


f 


1 


1 


^1 


10 


1 


^ 




400 



A R.A. (arcsec) 



Figure 12. From top to bottom; Column density map of C2H and 
C-C3H2, and their ratio, derived from the CARMA observations. 
For the column density derivation of both species, an excitation 
temperature of 20 K has been adopted. The other symbols are 
the same as in Figs. |8] and Fig. |9] 



core to core. Such values are much larger than the c osmic 
D/H ratio, estimated to be ~ 10"'' (Linsky et al. uOOq 'l. and 
comparable to those measured towards low-mass pre-stellar 
co res (R oueff et al. i2005i ) and infrared dark clouds (Pillai et 
al. [20071). This implies that the deuteration in the cores of 
120343 is as high as in colder and more quiescent environ- 
ments, despite the relatively higher gas temperature and tur- 
bulence, and confirms previous findings in other dense cores 



associated with intermedi ate- t o high-mass star- f orming re- 
gions (e.g. Busquet et al. I2OI0I . Pillai et al. I2OIII ). Interest- 
ingly, relatively high values are found also close to the UC 
Hll region associated with IRS 3, in MMB and IRS3-SW. 
Because NII2D is efficiently formed on dust grains, a strong 
UV radiation can heat the dust and cause NII2D evaporation 
thus increasing its abundance. On the other hand, a strong 
UV radiation could (at least partly) decrease the abundance 
of NII3, due to its interaction with H"*" to form NHj (e.g. 
Fuente et al. ll993l ). 

Concerning the NIl3-to-N2lI+ ratio, we find the largest 
value in MME (~ 80). The enhancement of the NH3-to- 
N2H"'" ratio can be understood when freeze-out of species 
he avier than He becomes important (see e.g. Flower et 
al. |200y), so that it is expected to increase when the star- 
less core gets closer to the onset of star formation. In this 
scenario, the fact that MME has the largest NH3-to-N2lI"'" 
ratio suggests that this core could be close to the onset of 
gravitational collapse, i.e. MME could be a cand idate mas- 
sive pre-stellar core. However, Palau et al. (|2007bl ) measured 
with the SMA a mass of only 0.7 M© from the 1.3 mm con- 
tinuum, while we find 23 M©. This discrepancy likely comes 
from extended fiux filtered out by the SMA, which means 
that the core is quite fiat and not centrally-peaked as ex- 
pected for a pre-stellar core. 

Based on the results of this work, we propose our fi- 
nal interpretation for the nature of each one of the 3 mm 
condensations: 

• MMA is probably a protostar candidate. Although it 
does not show any embedded infrared source, its relatively 
high Trot, large line broadening, and NH2D-to-NH3 lower 
than in other cores suggest that this condensation is evolved. 

• MME is likely a pre-stellar core, because it shows high 
NH2D-to-NH3 and NH3-to-N2H'*' ratios, is more quiescent 
than MMA and it does not appear fragmented into smaller 
condensatio ns when observed at higher angular resolution 
(Palau et al. l2007a ). Assuming a typical star formation effi- 
ciency of ~ 30%, the core, the mass of which is 23 Mq, has 
the potential to form an intermediate- to high-mass object. 

• the nature of MMB, MMC and MMD is less clear. Due 
to the low NH2D-to-NH3 and NH3-to-N2H+ ratios, MMD 
could be a protostellar object, consistent with clear hints of 
contraction motions seen in the NH3 (1,1) spectrum, while 
for MMB we found hints of expansion due to asymmetric 
emission in the two inner satellites (see Sect. I3.2.2p . Cer- 
tainly, all condensations are perturbed (MMB by the ion- 
isation front from IRS 3, MMC by IRS 1 and the outflow 
associated with it, MMD perhaps by a combination of both). 
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Only higher sensitivity and angular resolution observations 
will allow to better understand the nature of these cores. 



4.3 Interaction of IRS 1 and IRS 3 with the dense 
gas: an expanding cavity 

The most striking result of this work is the clear confirma- 
tion of a cavity opened by IRS 3 in the molecular surround- 
ing gas, and a tight interaction between this cavity and the 
surrounding dense gas. We have found several evidences of 
this: (i) the morphology of all the molecular tracers, espe- 
cially in C2H and C-C3H2, delineates a cavity around IRS 3 
and the 1.3 cm continuum map resolving the ionised gas 
perfectly matches the profile of the cavity; (ii) the NH3 in- 
tegrated intensity (2,2)/(l,l) ratios are large near IRS 3; (iii) 
the line widths are also large near IRS 3, specially in C2H 
and C-C3H2; (iv) the position- velocity plot of NH3 shows a 
U-structure typical of an expanding shell; (v) in the MMB 
core we found hints of expansion in the NH3 (1,1) spectrum 
due to different intensity of the two inner satellites. These 
evidences of such an interaction are shown for the first time 
in this work. 

If we put together all the results obtained, we specu- 
late about a possible scenario that describes the star forma- 
tion history in 120343: IRS 1 and IRS 3, both having bolo- 
metric luminosities of about 1000 Lq, seem to come from 
the same natal cloud while being clearly in different evo- 
lutionary stages, which points towards different generations 
of (intermediate- to high-mass) star formation in 120343. 
In this context, IRS 3 could have induced the formation of 
IRS 1, as could be inducing star formation on the west (in 
MM A). On the other hand, in this bright-rimmed cloud the 
star formation probably has not been triggered by the UV 
radiation from the Cygnus 0B2 association stars, because 
IRS 3, the massive star that formed first, is relatively distant 
from the bright rim, and the dense gas where we find the 
bulk of the current star formation activity is around IRS 3 
and away from the bright rim. Therefore, the star formation 
seems to be dominated by IRS 3, which has been caught in 
the act of pushing away and disrupting its natal cloud. 



5 SUMMARY AND CONCLUSIONS 

The protocluster associated with the centre of 120343 is 
an excellent location where the interaction between evolved 
intermediate- and high-mass protostellar objects and dense 
(starless) cores can be studied. We have derived the phys- 
ical and chemical properties of the dense gas by means of 
CARMA and VLA observations of the millimeter and cen- 
timeter continuum, and of several molecular tracers (C2II, c- 
C3II2, NII3, NII2D, N2II+). Below, we summarise the main 
results. 

• Morphologically, the dense gas is distributed in a fila- 
ment oriented east-west that passes in between IRS 1 and 
IRS 3, the two most massive and evolved objects. We resolve 
the dense gas into five millimeter continuum condensations. 
All of them show column densities consistent with poten- 
tially being the birthplace of intermediate- to high-mass ob- 
jects, although the masses derived from continuum suggest 
that they can form intermediate-mass objects more likely. 



• We confirm the presence of an expanding cavity driven 
by IRS 3 demonstrated mainly by the shape of the emission 
in the two FDR tracers C2II and C-C3II2, as well as by hints 
of expanding motions from both the position-velocity dia- 
grams and the asymmetric intensity of the two inner satel- 
lites of the NII3 (1,1) line of the millimetre core closest to 
IRS 3 (MMB). 

• The non-thermal line widths across the filament indi- 
cate that the gas kinematics is dominated by turbulence, 
similarly to other intermediate- and high-mass star-fornung 
regions and different from low-mass dense starless cores. 

• The rotation and kinetic temperatures derived from am- 
monia are on average larger than those typically found in 
cores associated with low-mass star-forming regions, espe- 
cially around the cavity walls. The most massive and ex- 
tended millimeter core, MME, shows physical and chemical 
signatures of an intermediate- to high-mass pre-stellar core 
candidate. 

• We have better constrained the spectral index of the 
radio-continuum emission associated with IRS 3, which 
turns out to be flat, and thus the ionised gas comes from 
a region photoionised by the B2 ZAMS star at the centre of 
IRS 3. 

• the column density ratio C2H/C-C3H2 is of the order of 
200-400 across the source and is higher where the dense gas 
is getting disrupted. 

• the deuterated fraction NH2D-to-NH3 is of the order 
of 0.1 in all cores, as large as in low-mass pre-stellar cores 
and infrared dark clouds. We find high levels of deuteration 
also close to the cavity driven by IRS 3. We speculate that 
evaporation of NH2D and NH3 destruction caused by the 
UV radiation field can influence this high deuteration. 

These findings undoubtedly confirm a tight interaction 
in 120343 between the most massive and evolved objects 
and the dense surrounding starless cores in several respects 
(kinematics, temperature, chemical gradients), and suggest 
that IRS 3 could have induced the formation of IRS 1 and of 
the other gaseous condensations accumulated on the cavity 
walls. However, the large-scale morphology of the molecu- 
lar tracers suggests that we are likely seeing only a limited 
portion of a very extended gaseous filament. Only a large 
pc-scale mosaic will allow us to fully trace the distribution 
of the dense gas in the region and delineate a complete view 
of the core population. 
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APPENDIX A: ANALYSIS FOR THE CASE OF 
DIFFERENT T^x FOR THE MAIN AND THE 

INNER SATELLITE LINES OF NH3 (1,1) 

Let us assume that the linewidth of the magnetic hyperfine 
components of the NH3 (1,1) inversion transition is large 
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Figure A-1. Geometry of MM A 



enough to make them unresolved, so that only the five elec- 
tric hyperfine lines, i.e. one main line ("m"), two inner satel- 
lites ("is"), and two outer satellites ("os"), are resolved. As- 
suming that the excitation temperature is Tex 3> Tbg, the 
ratio of intensities of the main line to the inner satellite line 



rMB(l,l;m) <rk(i-e- 



(A-4) 



giving that the optical depth of the main line must be Tm ^ 
3.1. Thus, taking into account that ria = 0.28 Tm, we obtain 



0.61 ^ 



1-e- 
resulting in 

Tcx(l, 



^1, 



0.23 ^ 



l;is) 



sC 0.38. 
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Tcx(l,l;m) 

The result is that the deeper layer traced by the satellite line 
is colder than the outer layer traced by the main line. If we 
assume that the kinetic temperature of MMA, obtained from 
the ratio Tm'b{'2, 2)/rMB(l, 1), is tracing the outer layer, the 
outer layer temperature is 22 K, while the inner layer tem- 
perature is between 5 K and 8 K. The higher temperature of 
the outer layer is indicative of external heating, as discussed 
in the text. 



TMB(l,l;is) _ Tex(l,l;is)(l-e- 



TMB(l,l;m) rex(l,l;m)(l-e— »)■ 



(A-1) 



The usual assumption is that we are observing a homo- 
geneous isothermal region, so that both excitation tem- 
peratures, rox(l,l;is) and Tcx(l,l;m), are equal, and that 
Tis = 0.28 Tm. In this case. 



(A-2) 



0.28 ^ TM^^y-^is) ^ 
rMB(l,l;m) " 

The lower limit corresponds to the optically thin case, while 
the upper limit is the optically thick case. 

In Sect. 3.2.2 we have shown NH3 (1, 1) spectra for each 
of the 3 mm continuum clumps. For the case of MMA, the 
intensity ratio of the inner satellites and the main line is 
0.23 ± 0.05, lower than the optically thin limit, 0.28. 

The assumption that the observed ratio is close to 0.28, 
and that the emission in MMA is optically thin, leads to 
inconsistent results. For an optical depth of the main line of 
Tm < 0.1, we obtain that Tex > 230 K, which is much higher 
than the kinetic temperature estimated from the intensities 
ratio rMB(2, 2)/Tmb(1, 1) (see text), Tl, = 22 K. This result 
is improbable, since we expect the excitation temperature 
to be, in general, lower than the kinetic temperature. 

The intensity of the (1, 1; m) is close to the kinetic tem- 
perature, indicating that the optical depth of the main line 
is probably Tm ^ 1. The optical depth of the satellite, how- 
ever, can be lower, so that both lines are tracing the emis- 
sion of different layers of the region observed: the main line, 
the outer layer facing the observer; and the satellite line, a 
deeper layer of material (see Fig. IA-1[I . The easiest explana- 
tion of the anomalous ratio rMB(l, l;is)/rMB(l, l;ni) is to 
assume that the region is not isothermal, and that the two 
layers at different physical depths, have different tempera- 
tures. So, the two excitation temperatures, rcx(l, l;is) and 
Tex(l, l;m), are not equal. Thus, 



Tex(l,l;is) 1 



Tox(l,l;m) 1-e— » 
Assuming that To; 



0.23. 
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rMB(l,l;m) 

(1, 1; m) is lower than Tk, 
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